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We consider the totally asymmetric simple exclusion process (TASEP) in discrete 
time with the sublattice parallel dynamics describing particles moving to the right 
on the one-dimensional infinite chain with equal hoping probabilities. Using sequen- 
tially two mappings, we show that the model is equivalent to the TASEP with the 
backward-ordered sequential update in the case when particles start and finish their 
motion not simultaneously. The Green functions are obtained exactly in a determi- 
nant form for different initial and final conditions. 
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I. INTRODUCTION 

The totally asymmetric simple exclusion process (TASEP) is a stochastic interacting 
particle system which serves as a paradigmatic model for nonequilibrium behaviour [1- 
3|]. The dynamics of this lattice gas model is characterized by the updating law. In one 
dimension on the integer lattice Z the most important cases of discrete-time updates are 
the backward-sequential, parallel and sublattice parallel updates For a finite number of 
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particles these dynamics can be defined through a master equation of the form 

P(x,t+1) = ^p,,,,P(x',t) 



(1) 



where x = {xj} describes the positions of the particles and Px.x' is the transition probability 
to go in one time step from a configuration x' to a configuration x. This transition probability 
is different for the various update schemes. For the backward-sequential update, each particle 
may take one step to the right with probability v if the target site is vacant at the beginning 
of the time step or becomes vacant at the end of the time step (due to motion of the 
particle in front). For the parallel update, the motion to the right is allowed only if the 
target site is vacant at the beginning of the time step. By iterating ([1]) one obtains the 
solution of the master equation for any given initial configuration x*^, i.e., the conditional 
probability to find a particle configuration x at time step t, given that the process started 
from configuration x°. This stochastic many-body dynamics have a natural interpretation 
in field-theoretic terms js, 5| where specific realizations of the process correspond to paths 
in the path integral representation of field theoretic quantities. Therefore, in analogy to 
the corresponding terminology in field theory, we refer to this time- dependent conditional 
transition probability as the Green function. 

For the first two cases, backward-sequential and parallel update, the Green functions of 



transition probabilities have been found 



systems defined on an infinite lattice 



yy explicit solution of the master equations for the 
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22|. The Green function has a determinantal 



representation similar to the one first discovered for the continuous-time definition of the 
process 17[ where particles jump independently after an exponentially distributed random 

-|3|]. This representation allows for a direct derivation of the current 



time with fixed rate 1 
distribution 



10| and has inspired a considerable amount of subsequent detailed analysis 

- Il4| and also of the 



of dynamical properties of the TASEP and related models, see e 
ASEP where particles are allowed to jump in both directions 
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16|. 
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type of discrete-time update, sublattice parallel, was first considered in some 



20|] and has subsequently been studied for various applications both analytically 



and numerically j^, |2l|. In this paper, we derive the Green function of the TASEP with 
sublattice parallel update which is defined as follows. 

Consider the process on Z, i.e. the one- dimensional infinite chain. Each site labeled 
by an integer i is occupied by at most one particle which can hop only to the right in a 
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FIG. 1: Space-time trajectories of four particles with appropriate weights on a chessboard. 

discrete time. At the first moment of time, we look at all {2i, 2i + 1) pairs. In each of them 
if the vertex 2i + 1 is free and the site 2i is occupied, the particle of the vertex 2i hops 
to the right with probability v and doesn't move with probability 1 — v. If both sites in a 
pair are occupied or empty or if site 2i is empty and site 2i + 1 occupied, the pair remains 
unchanged at that time step. At the next step of time we apply this rule of hopping to all 
pairs {2i + 1, 2z + 2). Continuing, we apply the updating rule to {2i, 2i + 1) pairs at each 
odd moment of time and to {2i + l,2i + 2) pairs at each even moment {23!. 



II. THE EQUIVALENCE OF THE TASEP WITH SUBLATTICE PARALLEL 
AND THE BACKWARD SEQUENTIAL UPDATES 

As noted above, the conditional probability to find particles at positions Xi < X2 < 
. . . < xn at discrete time t if these are in positions < < . . . < x% at initial moment 
of time is called the Green function of the process. The discrete space-time dynamics can 
be described by a set of trajectories on a triangle lattice which is obtained from the square 
lattice by adding a diagonal bond between the upper left and the lower right corners of 
each elementary square. Being occupied by an trajectory, diagonal bonds have a statistical 
weight V and vertical ones can have weights 1 or 1 — t;. It is convenient to draw trajectories 
of particles on a chessboard (Fig. [1]), where black rounds show initial positions of particles. 
We notice that diagonal bonds of trajectories can be located only on white squares. If we 
select a sublattice which contains upper left and lower right sites of white squares of the 
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FIG. 2: Space-time trajectories of particles starting and ending on the sublattice of white vertices 
(a) and rotated trajectories by 45 degrees clockwise (b). The time axis on the rotated lattice is 
directed down. 

chessboard denoted by white circles in Fig. [H we can see that particles effectively move on 
the sublattice of white vertices. There are some exceptions at the start and at the end of 
trajectories. Then, we have to consider four different cases to find a generalized determinant 
formula of the Green function. 

Consider first the case when space-time trajectories of particles start and end on 
the sublattice with white vertices, Fig. (the case of arbitrary initial conditions will be 
considered in the next section). If we choose initial points on the white vertices of the first 
row, coordinates of particles {x^} at initial times {T^ = 0} (i = 1, 2, . . . , A^) are even. 

The first transformation we use is a rotation of the set of trajectories by 7r/4 clockwise 
around the initial point {a; = 0, T = 0}. Considering the vertical axis as the new time 
coordinate, we obtain a new set of trajectories (Fig. [2b). This set represents a new discrete- 
time process on the square lattice with the unit time and space intervals corresponding to 
vertical and horizontal distances between neighboring sites. Starting points change their 
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FIG. 3: Final worldlines of particles after appropriate transformations, 
space-time coordinates as 

Tl' = x°/2, (2) 
xf = x\l2. 

Now vertical bonds have weights v and diagonal ones have weights 1 or 1 — f. We want 
to map them to the space-time paths of particles of the TASEP with backward-sequential 
update. To this end, we shift the coordinates in each row with respect to the previous 
above row by i — i- i + 1. Due to the second transformation, vertical and diagonal bonds are 
interchanged as it is shown in Fig. [31 The transformation of coordinates (in new units) can 
be written as 

(^",n = (^,^). (3) 

From Fig. [3] we see that worldlines on the transformed lattice represent trajectories of 
particles of the TASEP with backward-sequential update with initial space-time coordinates 
{x-"^, T/'°} and final coordinates {x-', T-"}, i = 1, 2, . . . , A^. The transition probability from 
space-time coordinates {x"°,T/"^} to {^■',T/'} is given by generalized determinant formula 



22| 



P(t"T"-t"T"- -t" t" It"" T"°- ^''O T"°- • t"° T"""! - Hpt A/f W (A^ 
where the matrix elements of A^ x A^ matrix M*^^) are 
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with the function {x,T) introduced in [22| 

FUx, T) = / dz{l - V + -f{l - zY'^z^-K (6) 

Substituting transformation ^ to the determinant formula (jlj) we obtain the Green function 
of the TASEP with sublattice parallel update 



P (xi, X2, . . . , xm\xI, X2, . . . , x^; r) = det M*-^-* 



(7) 



with matrix elements 



(TV) 



T + Xi ^ Xj^ 
2 



(8) 
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FIG. 4: Space-time trajectories of particles in the case, when the i-th worldline starts from the 
non-white vertex and ends on a white one (a) and rotated version of trajectories by 45 degrees 
clockwise (b). 



III. OTHER CASES OF STARTING AND ENDING POINTS 

Consider the case when the i-th particle starts its motion from an odd site and ends on 
an even one (Fig. H^). 
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As the weight of the first bond of the i-th trajectory is 1, we can set the beginning of 
motion at the nearest white (sublattice) site and then rotate the sublattice as in the previous 
case (Fig. Hb). Applying the shift transformation, we obtain for the initial coordinates of 
2-th particle: 



(a:?,T:°) = (x°,rx°/2l) (9) 

where \x] is the ceiling function. Substituting expressions of all starting and end points into 
the determinant formula (jl]), we derive the Green function for this case. 

The third case is when the i-th trajectory starts from the white (even) vertex and ends 
on a non-white one. We see, that if we add an additional vertical bond with weight 1 to the 
last node of that trajectory, the total weight of the whole path will not be changed. Then 
we can set the endpoint of i-th particle at time T + 1. Repeating two transformations, we 
obtain for coordinates of the end point: 

«,^") = (^.r^i) (10) 

The last case when the i-th trajectory starts and ends on non-white vertices is an obvious 
combination of the second and the third case. 

Generalizing all four cases of boundary conditions, we derive following transformations 
for the initial and final coordinates for all types of trajectories 

rr = rx°/2i, 
Ti' = r^i, 

< =Xi. (11) 

Substituting new coordinates (fTTj) into the determinant formula (jlj) we obtain the Green 
function of the TASEP with sublattice parallel update 

P{xi,X2,...,xn\x[,xI,...,x%;T) = det M^^\ (12) 

with matrix elements 

M^p = L - xi 1 - rf 1 ) • (13) 
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IV. DISCUSSION 

Having explicit determinant expressions for the Green function, we can compare their 
relative advantages and disadvantages for the three basic updates, the backward-sequential, 
parallel and sublattice-parallel one. Criteria for the comparison follow from practical use 
of the Green function in probabilistic calculations. To find a probability distribution for 
a selected particle or a correlation function for several particles in the TASEP, detailed 
information contained in function P (xi, X2, . • • , XTvIa^?, a;^, . . . , x^^r; T) should be reduced by 
summation over a part of the final coordinates {xi} for fixed initial coordinates {x^} (see 
e.g. Then, the first criterion for the comparison is simplicity of the summation 

procedure in different cases. The second criterion is simplicity of the matrix M^^'^ itself, 
because asymptotic calculations for large N and T need an elaborated analysis of resulting 



determinant expressions (see e.g. 13|, UM)- The third criterion is the presence or lack of 



particle-hole symmetry, which is essential for the derivation of single-particle probability 



distributions in some particular cases [lOj 



(A) The backward-sequential update. The form of the matrix elements M^^"^ in this case 



is especially simple 



Mf)=F,_,(x,-x;,r) (14) 

where function Fjn{x,T) is given by Eq.Q. The Green function 

P {xi,X2, . ■ . , Xiy\x^, X2, ■ ■ ■ , x%] T) is uniform in variables {xi}, so the summa- 



tion procedure is straightforward [lOj. A shortcoming of this update is lack of 
the particle-hole symmetry. Indeed, due to possible transitions for one time step 
Xj — )■ Xj + 1, Xj+i = Xj + 1 — )■ Xj+i + 1, . . . , Xj+fc = Xi^k-i + 1 — ^ Xi+k + 1, a hole can move in 



the opposite direction by jumps of length k > 1. 



(B) The parallel update. The form of the matrix M^^^ in this case is more complicated 
Q: 



=F,_,(x,-x°,r) (15) 

where 

F,.„(iV.T) = ± ± „ ^(.n + ;c + n-l). „ ^ 
^-^ ^-^ [k ^ n)\n\[m — ny. \ — v 

n=0 k=-n \ ' J \ I 
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The Green function obeys the particle-hole symmetry, but a drawback is in the determi- 
nant formula 

P{xi,X2,..., xn\xI 4,---, A; T) = {1- v)"" det M(^) (17) 

which depends on the number of pairs n of neighboring particles in the final configuration. 
Therefore, the sum over {xi} splits into groups by the number of clusters of connected 
particles. 

(C) The suhlattice parallel update. The Green function for this update is free of shortcom- 
ings of two previous cases. It is uniform in variables {xj}, obeys the particle-hole symmetry 
and has a relatively simple analytical form (??). A fee for this advantage is a rather com- 
plicated time dependence in (??) which involves both initial and final coordinates and the 
ceiling function \x'\ . Thus, we may conclude that a proper choice of the discrete time Green 
function strongly depends on peculiarities of the corresponding probabilistic problem. 
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